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ABSTRACT 



o 



, Traditionally, the degeneracy between the unknown radial profiles of total mass and velocity 

• anisotropy inherent in the spherical, stationary, non-streaming Jeans equation has been han- 

' died by assuming a mass profile and fitting models to the observed kinematical data. How- 

ever, mass profiles are still not well known: there are discrepancies in the inner slopes of the 
density profiles of halos found in dissipationless cosmological A^-body simulations, and the 
inclusion of gas alters significantly the inner slopes of both the total mass and the dark matter 
component. Here, the opposite approach is considered: the equation of anisotropic kinematic 
(-H ' projection is inverted for known arbitrary anisotropy to yield the space radial velocity disper- 

Q^l sion profile in terms of an integral involving the radial profiles of anisotropy and isotropic 

I ■ dynamical pressure (itself a single integral of observable quantities). Then, through the Jeans 

O ' equation, the mass profile of a spherical system is derived in terms of double integrals of 

^ , observable quantities. Single integral formulas for both deprojection and mass inversion are 

^ ■ provided for several simple anisotropy models (isotropic, radial, circular, general constant, 

Osipkov-Merritt, Mamon-Lokas and Diemand-Moore-Stadel). Tests of the mass inversion on 
NFW models with the first four of these anisotropy models yield accurate results in the case 
C*^ ' of perfect observational data, and typically better than 70% (in 4 cases out of 5) accurate 

^ I mass profiles for the sampling errors expected from current observational data on clusters of 

galaxies. For the NFW model with mildly increasing radial anisotropy, the mass is found to 
be insensitive to the adopted anisotropy profile at 7 scale radii and to the adopted anisotropy 
radius at 3 scale radii. This anisotropic mass inversion method is a useful complementary 
tool to analyze the mass and anisotropy profiles of spherical systems. It provides the practical 
means to lift the mass-anisotropy degeneracy in quasi-spherical systems such as globular clus- 
ters, round dwarf spheroidal and elliptical galaxies, as well as groups and clusters of galaxies, 
when the anisotropy of the tracer is expected to be linearly related to the slope of its density 
(Hansen & Moore 2006). 
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1 INTRODUCTION 

The major goal of kinematical modelling of a self-gravitating astrophysical system, observed at one instant, is to measure on one hand the 
total mass distribution (visible and dark matter), and on the other hand the three dimensional velocity streaming and dispersion moments. 
In other words, the modeller wishes to deduce the distributions of dark matter and of orbital shapes. The modeller has at his disposal, at 
best, maps of surface density (or surface brightness) and of the velocity field at each point, or else its moments (line-of-sight mean velocity, 
dispersion, skewness and kurtosis). 

The basi c equation for such kinematical modelling is the collisionless Boltzinann equation (hereafter CBE, but also often called Liouville 
or Vlasov, see lHenoi]! 19821) . which states the incompressibility of the system in 6-dimensional phase (position, velocity) space: 
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where $ is the gravitational potential (hereafter potential) and / is the distribution function, that is the density in phase space. Unfortunately, 
the resolution of the CBE is difficult, especially when projection equations are taken into accountQ In particular, the CBE presents a degen- 
eracy between the unknown potential and the unknown velocity field (given that observations usually limit the velocities to their projection 
along the line-of-sight, measured through redshifts). 

The traditional simpler approach has been to use the (first) velocity moments of the CBE, which are more easily related to observables, 
the Jeans equations that pertain to local dynamical equilibrium 

^ + (v- V)v = -V$-iv- (po-^) , (1) 

where p is the space density of the tracer used to observe the system, cr'^ is the tracer's dispersion tensor, whose elements are alj = 
Vi Vj — ViVj, and pcr^ is the anisotropic dynamical pressure tensor of the tracer. With the simplifying assumptions of stationarity and the 
absence of streaming motions, equation (TJ simplifies to the stationary non-streaming Jeans equations: 

V- (po-^) = -pV$ . (2) 

Using the stationary non-streaming Jeans equations (O, one can relate the orbital properties, contained in the pressure term with the mass 
distribution contained in the potential (through Poisson's equation). 

The small departures from circular symmetry of many astrophysical systems observed in projection, such as globular clusters, the 
rounder elliptical galaxies (classes EO to E2), and groups and clusters of galaxies, has encouraged dynamicists to assume spherical symmetry 
to perform the kinematical modelling. The stationary non-streaming spherical Jeans equation can then be simply written 

^jp^r) , „ /3 2 , . GMjr) 

___ + 2-pa, = -p(r)— 5— , (3) 
where M (r) is the total mass profile, while 

B(r) = 1 - '"^' + '^^ = 1 - ^ 

is the tracer's velocity anisotropy (hereafter, anisotropy) profile, with err = ^rr, etc., 00 = a^, by spherical symmetry, and with 13 = 1, 
0, —00 for radial, isotropic and circular orbits, respectively. The stationary non-streaming spherical Jeans equation provides an excellent 
estimate of t he mass profile, given all other 3D quantities, in slowly-evolving triaxial systems such as halos in dissipationless cosmological 
simulations 1 Tormen. Bouchet. & White! 19971) and elliptical galaxies formed by mergers of gas-rich spirals in dissipative A'^-body simulations 
jMamon et alj2006h . 

Again, one is left with having two unknown quantities, the radial profiles of mass and velocity anisotropy, linked by a single equation. 
In other words, we have to deal with a serious mass-anisotropy degeneracy. 

The simplest and most popular approach is to assume parametric forms for both the mass and anisotropy profiles. One can then express 
the product of the observable quantities: surface density profile E (fl) and line of sight squa re velocity dispersion profile af^g (i?) vs. projected 
radius R through the anisotropic kinematic projection equation l lBinnev & Mamonlll982h expressing the projected dynamical pressure P = 



PiR) = E{R)aL{R) = 2/ [(r'-i?') + 

J R 



"'^ (4) 



where equation lO is only valid for non-circular orbits, and where p = pa^ is the radial dynamical pressure. 

Inserting the radial pressur^ (eq. (5)) in the spherica l stationary Jeans equatio n ([3}, one determines the line of sight velocity disper- 
sions essentially through a double integration over p M dr. livtonon & LokasI l l2005bl. Appendix) have simplified the problem by writing the 
projected pressure as a single integral 

P{R) = E{R)aLiR) = 2G / J^p^oj [r, J?|/?(r)] pM — = 2 / K^^oAr, RWir)] P^c dr , (6) 



where they were able to determine simple analytical expressions for the dimensionless kernel ifpr oj for se veral popular analytical formula- 
tions of /3(r). Note that if one assumes isotropy, the equations are greatly simplified, and one finds l lTremain e et al. 1994; Prugniel & Simie^ 
fl997h ■K^proj (r, _R) = -v/ l — _R^/r^. Also, if /3 = est 7^ 0, the kernel can be expressed either in terms of incomple te Beta functions 
jMamon & Lokasll2005bl) . or in terms of the easier to compute regularized incomplete Beta functions jMamon & Lokasll200^) . With para- 
metric choices of the mass profile M{r) and anisotropy profile /3(r), one can fit for the free parameters of these two profiles that lead to the 
best match of the observed line of sight velocity dispersion profile. The drawback of this indirect method, even with the recent introduction 
of these simplifying kernels, is that the analysis is doubly-parametric, so that the derived parameters will be meaningless if one does not 
choose the correct form for both the mass and anisotropy profiles. 



Note that for dynamical studies, a fast computer code has been recently developed bv lAlard & Colombil lioo^ that solves the CBE in l-l-l dimensions, and 
a 3+3 dimension code is under development by G. Lavaux and these authors. 
^ For claiity, we hereafter drop the terai dynamical before pressure. 
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The next step in complexity is to perform a single -parametric analysis: either isotropy is assumed to directly obtain the mass profile, 
which we call the mass inversion, which is the focus of the present paper. Alt ernatively, a mass profile c an be assumed and one directly deter- 
mines the anisotro py profile through the anisotropy inversion, first derived bvlBinnev & MamonI (Il982[). with later and pro gressively simpler 
solutions found bv lTonrvl ( ll983h . lBicknell et alj ( ll989t) . ISolanes & Salvador-Solel ( Il990l) . and lOeionghe & MerritJ ( Il992h . One can attempt 
to lift the mass-anisotropy deg eneracy by cons idering together the variation with projected radius of the line-of-sight velo city di spersion and 
kurtosis fLokas..200Z ; .Lokas & Mamoij|200± For halos in cosmolos ical simulations, which are not far from spherical jjing & Suto, .20021 
and references therein) and nearly isotropic i Mamon & Lokasll2005bl an d references therein), viewed in p rojection, this dispersion-kurtosis 
analysis yields fairly accurate masses, concentrations and anisotropics jSanchis. Lokas. & Mamonll2004h . Unfortunately, the line-of-sight 
projection of the 4th order Jeans equation, required in the disp ersion-kurtosis me thod, is only possible when /3 = est, which does not appear 
to be realistic for elliptical galaxies formed by major mergers jPekel et alj2005l) . 

An even more sophisticated and general approach is to adopt a potential and minimize the residuals between the predicted and true 
observables, i.e. the distribution of objects in projected phase space (R,viob) (where iiios is the line of sight velocity) by one of several 
methods involving the distribution function: 

(i) A general global form for the distribution function is adopted, in terms of known integrals of motions. For example, in spherical 
systems with isotropic non-streaming velocities, the distribution function is a function of energy only, while in anisotropic non-streaming 
spherical systems it is a function of energy and the modulus of the angular momentum. Alas, t here is no known re alistic form for / = 
f{E, J) for anisotropic non-streaming spherical systems nor for non-spherical systems, although Iwoitak et al.l ( l2008h have recently shown 
that cosmological halos have distribution functions that can be written f{E, J) = j2{/3~-/3o) j_ j'^/j^^ \ where we adopt 
hereafter the notations /3o = /3(0) and /3oo = lim, — ,00 /3, where Jo is a free parameter related to the 'anisotropy' radius where /3(r) = 
(/9o + /3oo)/2. Unfortunatelv. lWoitak et al] do not provide an analytical formula for fsiE). 

(ii) A set of elementary distribution functions of E ov {E, J) is chosen, as first proposed bv lPeionghd ( Il989l ). then lMerritt & Sah3 ( Il993h . 
and applied to elliptical galaxies by Gerhard et al. ( 1998). One then searches the linear combination of these distribution functions, with 
positive weights (to ensure a positive global distribution function) that minimizes the residuals between the predicted and true observables. 
However, there is no guarantee that the set of elementary distribution functions constitute a basis set, so that some global realistic distribution 
functions may be missed. Moreover, the distribution function may depend on an additional unknown integral of motion. 

(iii) A set of delta- distribution functi ons, / = f(E,J) i s chosen, in other words one considers orbits of given E and J jSchwarzschildl 
11979 : Richstone & Tremaine 1984; Sve r&Tremainelll996h . A gain one searches for a linear combination of these orbits that minimises 
the residuals between predict ed and true ob servables, again enforcing positive weights. These weights are obtained eith er by averaging 
the observables over an orbit jSchwarzschild) or by continuously updating them jSver & Tremainel;lde Lorenzi et al.ll2C)07[). This method i s 
powerful enough to handle non-spherical potentials. Despite concerns about convergence jCretton & Emsellenill2004 IValluri et al.ll2004h . 



the orbit-superposition method, if properly implemented, does reproduce the correct solutions ( iRichstone et alJl2004l ; lThomas et al.ll2004l) 

The potential can be adapted from the observations, assuming constant mass-to-light ratio (M/L) if the observed density is a surface 
brightness, or constant mass-to-number ratio (M/N) if the observed density is a surface number density. If spherical symmetry is assumed, 
this involves a choice of M/L or M/N, the deprojection of the surface density map, and then Poisson's equation is easily inverted to obtain 
the potential from the density. For axisymmetric syste ms, one can deproject the surface density maps into a potential assuming it to be the 
sum of g aussians l[E msellem, Monnet. & Baconlll994h . One can add to the potential a possible dark component given in parametric form 
(see, e.g.. lWilliains et aL2009i) . 

Alternatively, instead of using distribution functions, one can fit the distribution of objects in projected phase space by the multiple 
parametric a djustment of the mass and anisotropy profiles, as well as pos sibly the velocity distribution in space (which could be non- 
gaussian, see lKazantzidis et alj|2004l : IWojtak et al]|200^ : [Hansen et alj|2006h . as in the MAMPOSSt method (Mamon, Biviano & Boue, in 
preparation). 

Returning to direct single-parametric estimatio ns, the mass profile of astronomi cal systems does not seem to be better established than 
the anisotropy profile. Indeed, despite early claims l lNavarro, Frenk. & Whitjl 19961 hereafter NFW) of a universal density profile for the 
structures {halos) in dissipationless cosmological A'^-body simulations of a flat Universe of cold dark matter with a cosmological constant 
(hereafter A CDM), there has been an ongoing debate on whether the inner slope is steeper (jpukushige & Makino 1997; Moore et al. 199^ 
or shallower jNavarro et al]|2004IStoehr et al.ll2002l;IStoehill2006l ). Furthermore, the inclusion of gas in cosmological simulations can lead 
to much steeper dark matter density profiles "l lGnedin et alj|2004l) . Indeed, the dissipative nature of baryons leads them to accumulate in the 
centers of systems, not only in spiral galaxies, as is well known, but also in elliptical galaxies, for otherwise the NFW-like mass distributi on as 
found in ACDM halos would lead to a local M/L and aperture velocity dispersion much lower than observed ("m amon & Lokas' 2005a|), and 
the do minance of baryons in the center and dark matter in the envelopes has been recently confirmed by X-ray measurements ( Humphr ey et al.l 
l2006h . Moreover, the dark matter dynamically re sponds to the baryons that dominate in the i nner regions, to reach steeper slopes than they 
would have had without the presence of baryons ( lBlumenthaletal.ll 19861 ; lOnedin et al]|2004l) . But the final density profile of dark matter is 
expected to be very sensitive to the details of the baryonic feedback processes. 

On the other han d, the a nisotropy profiles of the halos in dissipatio nless cosmological simulations appears to be fairly universal (see 
the compilation by iMamon & L okas 2005 b) and references therein, andTwojtak et al.ll2008l) . although galaxy-mass halos have somewhat 
more radial orbits than cluster-sized halos lAscasibar & Gotfl6bej|2008l) . Also a similar shape of anisotropy profile holds in A'^ -body-l-SPH 
simul ations of merging spirals galaxies, including gas, but with a ratio of anisotropy to virial radius that is ten times smaller jDekel et"al] 
I2OO5I) . Moreover, dissipati onless A^-body sim ulations (cosmological and binary mergers) indicate that the anisotropy is linearly related to the 
slope of the density profile l lHansen & Moore 2006). although the trend is less clear in elliptical galaxies formed in A^-body-l-SPH simulations 
as merger remnants of spiral galaxies (Mamon et al. .2006 ). because of the dynamical interaction of the stellar, dark matter and dissipative 
gas components. 

In this paper, we derive and test the mathematics of the mass inversion. We begin in Sect. 12. II with a reminder on the kinematic 
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deprojection of isotropic systems, followed by the mass inversion of isotropic systems in Sect. 12.21 We then develop in Sect. l2.3l our algorithm 
for the kinematic deprojection of anisotropic systems, and in Sect. 12.41 we deduce the mass profile with the Jeans equation ([3}. In Sect.[3] we 
test our mass inversion methods. 

The reader in a hurry might want to skip the mathematical details. (S)he will find the general anisotropic deprojection formulae in 
equation l |32| l. with special cases given in equations ( |39| l [radial orbits], l l44b [circular orbits], and in equations ( |62| l [constant /3 < 1, 
Osipkov-Merritt, Mamon-Lokas, and Diemand-Moore-Stadel], with Ca given in Table [2] and kernels Ki given in equations ([37} [constant 
/3 < 1], l l48b [Osipkov-Merritt], l l53l l [Mamon-Lokas], and ( I61t [Diemand-Moore-Stadel]. The formulae for the mass inversion will be found 
in equations ( 1691 ) [general], MA\ [radial], ( |76| ( [circular], and J83b [constant /3 < 1, Osipkov-Merritt, Mamon-Lokas, and Diemand-Moore- 
Stadel] with the same and kernels, and with Dp also given in Table [2l 

In the very late stages of this work, w e came acro ss a draft o f lWolf etaljfeoogi) , who independently developed an analogous method for 
anisotropic kinematic deprojection. While IWolf et al.l produce a general formula for kinematic deprojection, the present article also provides 
simpler formulae for the kinematic deprojection with specific simple anisotropy profiles, as well as general and specific formulae for the 
mass profile. 



E(i?)= / pWd. = 2 / 7f7^> (7) 



2 METHOD 

2.1 Kinematic deprojection of isotropic systems 

We begin by reviewing the mathematical formalism for the kinematic deprojection of isotropic systems. The structural projection equation, 
relating the space density p(r) to the (projected) surface density E(_R); 

p(r) r dr 

is inverted through the usual Abel transform, whose derivation we recall in appendix [A] as we will use it in the following subsection. One 
then recovers the well-known structural deprojection or Abel inversion equation 

1 r AH AR 

In the case of isotropic velocities one can express the projected dynamical pressure E crf^j, in terms of the dynamical pressure p with 
the isotropic kinematical projection equation, obtained by setting /3 = in the anisotropic kinematic projection equation (|5j: 

2 rAr 

, . (9) 

Equation (|9) is the strict analog to equation (|7]l, where the tracer densityp is replaced by the dynamical pressure p — pa^ and the surface 
density E is replaced by projected pressure P — T} af^^^With these replacements, the structural deprojection equation ([8j turns into the 
isotropic kinematical deprojection equation 

pUr)^[pir)a^ir)W_.^-- ^ ^_ ■ (10) 



2.2 Mass inversion of isotropic systems 

Now, from the stationary non-streaming spherical Jeans equation ((S}, with the isotropic condition ((5 = 0), the total mass profile is trivially 

^' Gp Ar TvGpArJ^ AR ^R2 _ ^ ' 

With the variable substitution R = r u, we can avoid the singularity in the surface term of the derivative of the integral of equation dlOt or 
dl It by writing 

' ( ^ = ^ = _iA r P"( ) '^^'^ ^ r ^^iAR_ 

Ar nArJ^. AR ^rT^ nj^ Vil^ 7vr di?^ y^2^ ' ('^> 

where P"{R) = d^P/AR^. Inserting the right-hand- side of equation ( 112b into the first equality of equation dl lb . we then obtain the isotropic 
mass inversion equation 

A^ (EaL) RdR 



r 



, , AR^ V-R2 _ r2 

M{r) = — — , (13) 

^ ' G f °° AE AR 



AR - r2 



Given the isotropy, the space velocity dispersion is equal to the radial velocity dispersion, so we drop the subscript V. 
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Table 1. Nomenclature 



Definition 


full expression 


abbreviated expression 


space radius 




r 


projected radius 




R 


projected pressure 


los 


P 


radial pressure 


po-2 


P 


tangential pressure 




pg = 


anisotropy 






circular velocity 


x/GM(r)/r 





where we used the structural deprojection equation ([8j to replace the density in the denominator. The isotropic mass inversion equation can 
be further simplified, expressing the circular velocity, v1 = GM/r afl 



vl{r) = 



Tvp{r) 



d^P RdR 



(14) 



Unfortunately, the mass and circular velocity profiles require the second derivative of the (observable) projected pressure P — 'Saf^glR). The 
singularity (^R^ — r^^ in the numerators of equations il3l and l ll4t prevents one from expressing the mass profile with single integrals 
derivatives of the projected dynamical pressure after a suitable integration by parts. 



2.3 Kinematic deprojection of systems of arbitrary known anisotropy profile 

2.3.1 General anisotropy 

The anisotropic kinematic projection equation (O is strictly valid for non-circular orbits (finite /3). For circular orbits (or = 0) equation ^ 
yields 



where 

Pe = (1 - P)P = pc^e 

is the tangential dynamical pressure. To guide the reader, Table[T]reviews the nomenclature adopted in this paper. 
We repeat the steps of the standard (isotropic) Abel inversion (Appendix 0, now defining 



J(r-) = 



RdR 

- 7-2 



dP 
dR 



^/W^^dR, 



where equation J17t is obtained by integration by parts (the surface term is for P{R) cc R " with a > 1). 

For non-circular orbits, inserting the projected pressure (eq. (5)) into the definition of J (eq. 1161 ). one finds 



J(r) = 2 

2 / psds 



RdR 

RdR 



s ds 



r2 \ 1 

+ 1 d psds . 



ds 



R^ dR 



(15) 



(16) 
(17) 

(18) 
(19) 
(20) 



where equation l |19t is obtained after reversing the order of integration and the two inner integrals of eq. ( |19t are worth 7r/2 and (7r/4)(r2 + 
s^), respectively. Differentiating equation J20b . one has 



dJ 
dr 



/ I3p — 



(21) 



^ We've never encountered in the literature the mass profile written in this direct fashion, although lRomanowskv et al] i2009l) gave the equivalent expression 
-2 Z""" d{P'/R) dR 



Kir) 



TTp{r) 



dR y/W 



The expression in equation I14t seems preferable as the differentiation is performed in a single pass. 
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Now, equation ( 117b can be differentiated to yield 

dj r dP dR 

where the second equality in equation ( I22t comes from equation l llOt . Equations l l21t and V.2\ yield 

/■°° ds 

w(r)=Pi.o(r)- / /3p— . (23) 



Equation ( I23t is an implicit integral equation for p with piso (eq. 1101 ) and (3 known. For finite /3 < 1, we solve for p by differentiating 
equation ( I23t . to get the differential equation 

rl^p^jL^ (24) 
1-/3 r 1-/3 



Now, if we write 

, r0 +13 p _1 d[gp) 
1 — /3 r (; dr ' 

then equations l l24t and l |25t lead to 



(25) 



1 / ' 

where the upper limit at infinity ensures that the radial pressure p — pa^ does not reach negative values at a finite radial distance. But 
equation l l25t directly gives 

din/ _ -r/3'+/3 
dlnr ^ 1-/3 ' 

hence 

9(r) = g{ri) exp - / — (27) 



1-/3 s 

for any arbitrary ri. With equation J27t . equation l l26t leads to 



1-/5 s ; - V y,,^ 1-/3 * ; 1 - /3 

/ f"ti3' + f3dt\ pL„ , 

where the second equality is obtained adopting ri — r. 

One wishes to avoid the second derivative of the observables that occurs in the expression of equation ( 112b for p'is^oiT), which will 
amplify any uncertainties on the measurements of these observables. Integrating by parts the integral in equation ( |28b . we finally obtain 

where piso is given in equation dlOb . and where 

. , X / r tp' +13 dt\ 



which is provided in Table|2]for various simple anisotropy models. 

One may prefer to use the tangential dynamical pressure instead of the radial one, as it can be expressed in a slightly simpler form: 

Peir) = exp^-^ Y^T^^''""'^^' ^^^^ 

, , n /3 f f 13 dt\ ds 

= Piso(r-) - / Piso:; ^ exp ' ' ' 



Piso(s)B/3(r,s)^ , (32) 
as similarly derived in Appendix IB] and where 
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Table 2. Terms in equations 134) , (62) and (67) for specific anisotropy profiles 



Anisotropy model 


-4,3 (r, s) 


B;3(r,s) 


C,fl(r) 


Dp{r) 




(eq. (D) 


(eq. (33)) 


(eq. (62)) 


(eq. (13) 



^ = est 



Osipkov-Merritt (eq. 1451 ) 



Mamon-Lokas (eq. 1491 ) 



Diemand-Moore-Stadel (eq. (58)) 



p /^x,3/{i-/3) p ^ (3-2/3)/3 



/3 /^N/3/ii-pj ^ 

1-/3 UJ UJ 



2 1-/3 1-/3 

.,2 /„2_„2\ /„\2 /„2_„2\ ^ /„\ 2 ^2 _|_ 5q^2 



\a) + a2 '^'^P y 2a^ J \a) "^^^ \~2a?~ j a \a) 

r + a/ s\ (r + 2a)s /r\r + 2a 2r 

s + a Vs + 2a/ (s + 2a)2 Va/ a a + r 



- a" 



,1/3 A Z_ .,1/3 J 



2 /r-\ 1/3 Sai/^* - Sr^/^* 



(^ai/3 - ri/s)" (ai/3-ri/3)^ (ai/3-ri/3)^ sVa/ ai/3 - ri/3 



Notes: the Diemand-Moore-Stadel values are restricted to r < a. 

which is provided again in Table|2]for our simple anisotropy models. The radial pressure is then simply p(r) = pB{r)/[l — Pij-)]. 

The expressions for the dynamical pressure (radial or tangential) are made of single integrals involving piso, which is a single integral 
itself. Hence, the dynamical pressure is expressed in terms of double integrals. For simple anisotropy profiles, we can simplify the dynamical 
pressure to single integrals by inserting the expression forpiso(s) (eq. 1101 ) in equation ( I32t and inverting the order of integration. This yields 

1 P dP /? f f p dt\ ds 

Pe{r) = piso(r) + - / X5 di? / ^ exp ' 



7V dR 7,. 1-/3 "V 7, 1-/3 t I sVn^ 



= Piso(r) + - J P'{R)dR j Bp{r,s)—^^= (34) 



and for simple /3(r ), the inner integral can be expressed in closed form, as we shall now see. 



2.3.2 Case of finite (3 = est < 1 

Equation i29\ with Ap from Table [2] leads to 

Piso(r) /3 



p{r) 



,/3/(l-/3) 



l-/3(r) (1-/3)2 



Piso S 



-/3/(l-/3) 



ds 



(35) 



Using equation 04l ( with B/s from Table|2l one obtains a single integral representation for the tangential dynamical pressure 



/ ^ ^ \ , ^ P /3/(i-/3) 
pe(r) = piso(r)H z -r 



TT 1 



Piso(r-) + — 



-(3' 
P 



dP 
dR 



dR 



-/3/(l-/3) 



ds 



1 

2tv 1 - /? r 

where the second equality of equation ( 
is 



^ A'cst {^] dR 



( f 



dR 



\R 



(36) 



1 

u , -, — 



is obtained with the change of variable t — 1 — s^/i?^. The dimensionless kernel in equation ([36] 

/3/2 



(37) 



2' 1 

where B{x, a, b) = t°~i (1 — t)''^idf is the incomplete Beta function. Integrating by parts the integral in equation l l36t , we finally obtain 
after some algebra a single integral expression for the tangential pressure that does not depend on derivatives of the observations: 



Pe{r) =Piso{r) + 



2^ {1-pf 



P{R)R 



-{2-/3)/ {1-13) 



B 1 



R^'2'~l-/3 



dii- 2(1-/3) 



P{R) dR 



, - r-2 R 



(38) 



The surface term R ^^^^ B [l - r^/R^, 1/2, -/3/(2 (1-/3))] P(i?) that occurs in the integration by parts goes to as i? ^ cxd. Indeed, 
for X = r/R and c = — /3/2/(l — /3), one has x^~^'^ B(l — x^, 1/2, c) = —x/c + 0{x'^) and moreover P — » 0. In practice, if a 
progranmiing language does not provide the incomplete Beta function, but only the regularized incomplete Beta function, as I{x, a, b) = 
B{x, a, b)/B{a, b) = T{a + b)B{x, a, b) /[r{a)T{b)], one should then be careful that T{b) diverges when the last term b in the incomplete 
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Beta function is a negative integer, i.e. when (3 — 2n/ (2n + 1) = 2/3, 4/5, 6/7... (n being a positive integer). Luckily, B{x, a, b) always 
converges to finite values [f| 

In the limit P —> everywhere, equations ( 1351 ), ( 1361 ), and ( |38t all reduce to p{r) — Piso(r), as expected. 

2.3.3 Case of radial orbits: (3 — 1 

For radial orbits, differentiation of equation J23t leads to 

pir) = -rpUr) = - I ^ • (39) 

2.3.4 Case of circular orbits: (3 — cxj 

For circular orbits, we proceed in a similar fashion: inserting the projected pressure (eq. 1151 ) into the definition of J (eq. 1161 ), one finds 



/•oo i 

" ' V+.^) ^i^, (40) 



2 
and 

— = -7rr I p9 - / p9 — 1 . (41) 

Equations ( I22t and ( I41t lead to 

d / 1 dJ\ , pe , 1 d(/pe) 

-^T- — -T- =PeH = Piso = 7 — 1 , (42) 

dr \Tvr dr / r j dr 

whose solution is given by / = r: 

where the 2nd equality is found by integration by parts, for which the surface term, limr^oo ''Piso(''), vanishes for dlnp/dlnr + 
din Af/dlnr < (as derived from the Jeans equation |[3l), as is the case for reasonable mass and tracer density profiles. Inserting piso 
(eq. mOI ) into equation ( 143 1 and inverting the order of integration, we finally obtain the single integral expression for the tangential pressure: 



1 Z""" dP r 1 1 -1 / r \ 

po(r) = / — -r- — , H — cos — 

' TV dR ^R2 _ ^2 r \rJ 



dR . (44) 



2.3.5 Case of Osipkov-Merritt anisotropy 

For the Osipkov-Merritt ( IOsipkov|[T97^ : lMerritll 198^ anisotropy 



2 



= ^ . (45) 

^2 _|_ ^2 

equations ( 13 It reduces to 

Peir) = - ^ exp (~ \''J ^ pLods = (pcr,').^^ ('^) " ^ ^ exp (^^-^^) Pisosds , (46) 

where the last equality is again obtained after integration by parts or from equation ( I32t . Equation ([34} yields (see Table|2} a single integral 
representation for the tangential dynamical pressure: 

/ N 1 f r'\ r <^P .r, f" f sds 

peir) = Piso(r) + — r exp — ^ / —dR 



7ra2 V2a2 / / dR " V 2a 





where the dimensionless kernel Kqm is 

A'oM(?i) = V2F{u) = exp (-u^) erfiw , (48) 
^ An SM macro for B(x, a, b) is available upon request. 
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where 



F(n) 



exp (— u^) erfii 



^/^T exp (—u^) erf(m) 



2i 



is Dawson's integral and where erfi(a;) is the imaginary error function. Note that Dawson's integral is available in most software packages for 
mathematical functions [f| Equation ( 147b can also be found by inserting the expression forpiso(s) (eq. 1101 ) into equation i46i and reversing 
the order of integration. 

2.3.6 Case of Mamon-Lokas anisotropy 

For the simple anisotropy profile that lMamon & LokasI ( l2005bl) found to fit well ACDM halos 
1 r 



one obtains 

pe{r) = -{r + 2a) 



ds 



ISO I r» 

s + 2a 



= Piso(r) -{r + 2a) 



ds 



"° (s + 2a)2 



(49) 



(50) 



where the first equality is from equation ( 131b . while the second one is obtained after integration by parts or from equation ( 132b . Equation ( 134b 
now yields (see Table|2j the single integral expression for the tangential dynamical pressure: 

, , r/2 + a , , 2, , , , 1 , „ . dP da 



1 r + 2a 
Piso(r) H 



dn \ a a J 



TV a" 

where the dimensionless kernel Kuh, using X — R/a, x = r/a and y — s/a, is 



(51) 



KMLiX,x) 



dy 



1 

X2 



y/X^-y^ (3/ + 2)2 



de 



{cos9 + 2/Xy 



1 VX^ - a;2 



4-^2 2 + x 
_i_ {4 + x)V2^ 

12 (2 + 1)3/2 



VX2 



X2 



2 + X 



(4~X2) 



(X2 - if/^ 



3/2 



tan 



2-X X-x 



2 + X y x + x 



tanh 



X-2 X~x 



x + 2 y X + x 



for X < 2 , 
for X = 2 , 
for X > 2 . 



(52) 



(53) 



where equation ( |52b is found through the variable substitution y — X cos 9. Equations (jSTJ and l |52b can also be found by inserting the 
expression for piso (eq. 1101 ) into equation ( 150b and reversing the order of integration. 



2.3. 7 Case of generalized Mamon-Lokas anisotropy 

The velocity anisotropics in halos in cosmological A*' body simulations do not always fit th e Mamon-Lokas formula (eq. II49I ). but instead, 
j3{r) shows h alo to halo variat ions in its limits at r = and r oo dWoitak et al1l2008l) . Hence, a more general form for the anisotropy 
profile is (e.g. lTiret et al.ll2007l) 



/3(r) =/3o + (/3oo-/3o) 



r + a 



(54) 



The Mamon-Lokas anisotropy is the special case with /3o — and /3oo = 1/2. For /3o < 1 and /3oo < 1, inserting equation ( 154b into 
equation ( 131b yields, after some algebra: 

,-/3o/(i-^o) [(1 _ , + (1 _ p^) ^fo/(i-M-0^/a-0^) pj^^ . (55) 
For /3o < Poo = 1, the same procedure gives 



SM macros for crfi(x') and Dawson's F{x) are available upon request. 
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I \ ( \ f*o/(l-)3o) 

P0{r) = — exp r*^"'^ ™' 

\ 1 - PO / 

For /3oo < /3o = 1 III we similarly obtain 



P6)(r) = -exp ( I ^' 



exp 



exp 



s/a 
a/s 



-/3o/{l~/3o) ' 



-/3o„/(l-/3oo) / 



Piso ds . 



(56) 



(57) 



The integrals of equations ( 155b . J56b . and l l57t . are essentially double integrals, since they involve pJ^Q (eq. 1121 ). Single integral solutions do 
not appear to be possible to derive, even for the simple case of /3o = (unless j3oo ~ 1/2, i.e. the Mamon-Lokas anisotropy model). 



2.3.8 Case of Diemand-Moore-Stadel anisotropy 

Finally for the other simple anisotropy profile that lDiemand. Moore. & Stadell ( |2004| . Sect. 3.3.2) also found to fit well ACDM halos 



/3(r) 



1/3 



r < a , 
1 r > a , 



(58) 



we obtain 

pe{r) =piso(r-)- 

p{r) =Prad(r) 



(ai/3 _ri/3)'^ 



/ 1/3 1/3 



2 ds 

' 1^ 



for r < a . 
foi r > a . 



(59) 



where equation i59\ is obtained from equation l l32t . while the equation ( |59b comes from the pure radial orbits for r ^ a (eq. 1581 ). Again, 
for r < a, the integral in equation l |59t is essentially a double integral (because of p,'ao)i and a single integral solution can be obtained using 
equation ( |34t . yielding (with Table|2} 



Peir) = PisoW + 



Piso(rj 



1 



1 



dP ,^ r («^'^ - ^''') 



/3 - S" 



:ds 



where the dimensionless kernel is 



s 

_R r 



dP 

-t^Kdms (— , - ) dR, 
aft \ a ay 



(60) 



(■cos"l(a;/X) 



A2/3 





.£(1/6) 
r(2/3) 

for A — R/a, x = r/a and y = s/a. 



s2/3g Al/3 



-l(x/X) 



AO 



COsl/3 



+ 



:os~^(a;/X) 



B 



(- - -)] 


A-2/3 


\X'^' 6' 2 J 


2 



^ r(i/3) f 1 1 

r(5/6) Ia2'3'2 



de 

^-1/3 



_1 / 1 



a; ' 



(61) 



2.3.9 General expression for the tangential pressure for specific anisotropy profiles 

The expressions for the tangential pressure for the cases of constant, Osipkov-Merritt, Mamon-Lokas, and Diemand-Moore-Stadel anisotropy 
(eqs. 1361 , |47| , 1511 , and 1601 , respectively) can all be written in the form 



, ^ 1 ^ / N /■ dP ^ „ 1 r dP 



dP, 



(62) 



where the second equality of equation ([62} is found with equation dlOb and where Cfl(r) and KpiX^x) are dimensionless functions such 
that 



\a a. 



Bi3{r,s) dy 



Vr^ — y 



(63) 



with Ci3{r) given in Table [21 and given in equations {37}, ([48}, ([53} and ( I61t , respectively. For the Diemand-Moore-Stadel anisotropy 
model, the upper integration limits in equation i62\ should be replaced by the anisotropy radius a. The second equality of equation ( 162b 
allows the kinematic deprojection with a unique single integral. 



Decreasing anisotropy profiles are found for some regular halos ^Woitak et al1l2008h . although the central anisotropy is never unity. 
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2.4 Mass profiles of spherical systems with arbitrary known anisotropy 

2.4.1 General mass profile 

The mass profile is obtained through stationary non-streaming spherical Jeans equation (O, which writes 



GM , 2 



(64) 



Now, equation ( 124b reads 

P 



pLo P' + 13 p _ pi^, 



+ 



rP' + l3 pe 



1-13 1-/3 r 1-/3 (i_/3)2 r ■ 
Inserting p' from equation ( 165b into equation ( 164b yields the general mass inversion equation (dropping the dependencies on r for clarity): 



(65) 



-(l-/3)p^(r)=pL(r) + 

where the dimensionless function 

rd/3/dr + (3-2/3)/3 



l3' + {3~2f3) (3/r 
1-13 



Pe(.r) ^Pisa(r) H pe{r) , 



D0{r) = 



1-13 



(66) 



(67) 



is given in Table |2] for four anisotropy models. Inserting the general expression for pe into equation ( 166b . and converting the mass into the 
circular velocity with v^{r) = G AI{r) jr gives either 



[1 - /3(r)] p(r) vi(r) = -rp-^^^{r) + 
(from eq. 1291 ) or 



r/3' + (3-2/3)/3' 



■Si 



Piso 



\l-P(r)\p{T)vl{r) 



r/3' + (3-2/3)/3 



1-/3 



exp 



1-/3 

/3 

-/3 t 



exp 



/3 dt\ ds 

s 



1-/3 f 



' Piso 



(68) 



(69) 



(from eq. 1281 ). Alas, both forms (eqs. 1681 and 1691 ) involve the second derivative of the observable P, hence the second form (eq. 1691 ) 
seems preferable to use as it is simpler. However, for simple anisotropy profiles, the double integral of equations ( 168b and ( 169b can be 
simplified to single integrals, or equivalently, single integral expressions for pe exist, which can be inserted into equation ( 166b to obtain a 
single integral expression for the mass profile. 



2.4.2 Case of isotropic systems 

For isotropic systems (/3 = 0), equation ( 



I trivially leads to 



p{r) 



which is equivalent to the first equality of equation dl lb . 



(70) 



2.4.3 Case of finite (3 — est < 1 

For finite (3 — est < 1, while equation ( |69b becomes (with Dp from Table[2]l 



GM 

(1-/3) P^^(r-) =Pis, 



/3(3-2/3)/r 
1-/3 



^/3/(l-/3) 



Piso ds . 



a single integral expression is found inserting the tangential pressure (eq. 1381 ) into equation ( 166b to yield 



(1-/3) p—(r-) = pL + 



/3(3-2/3)/r 



1-/3 



Piso + 



/3 



27r (1-/3)2 



p(^)^-(^«/(i-«^/^_rJ 1 di?- 2(1-/3)/'°" 



i?2'2' 1-/3 



/R2T772 



(71) 



.(72) 



2.4.4 Case of radial orbits: 13 — 1 

For radial anisotropy, equations ( |39b and ( 164b simply yield 

GM(r) 



p(r) 



= 3pieo + rpiso 



However, using the change of variables R = r cosh it, the last equality of equation J39b yields 
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I Prad , 1 
Prad = \ 

r 7rr 



dR 



hence, from equation l |64| l 
1 



TV p[r) 



3P +RP 



RdR 



(73) 



(74) 



Alas, equation f!4\ involves a triple differentiation of the observables. 



2.4.5 Case of circular orbits: (5 —oo 

For circular orbits, the first term in the stationary non-streaming spherical Jeans equation lO vanishes, and one is left with the trivial relation 



p{r) 



Piso + 



Piso ds 



(75) 



where we made use of equation l l43t for the last equality. Integrating the last equation by parts, or equivalently, using equation we get 

1 



2/ ^ 2 



TV p{r) 



dP 
dR 



H — cos 



R 



dR. 



(76) 



2.4.6 Case of Osipkov-Merritt anisotropy 

For Osipkov-Merritt anisotropy (eq. 1451 ), equation l |66t leads to (with D/s from TableO 



-(r) = 



- Piso + 



+ 4 r 



exp 



Piso ds . 



\ / ^ r \ 

A single integral solution to the mass profile is obtained by inserting pg of equation l l47t into equation J66b , yielding 

R^ 



GM , , 



, pLo + {r^ + 5a^) 



1 / 



dP 
'dR 



exp 



(77) 



(78) 



2.4. 7 Case of Mamon-tokas anisotropy 

For lMamon & LokasI anisotropy (eq. 1491 ), equation l l66t brings (with Dp from Table|2j 

, ds 



GM r, r + a , 

P^-ir) = -2 Piso + 4 

r + 2a 



PiE 



s + 2a 

The single integral solution, found by inserting pg from equation iSll into equation l l66b . is 

GM, ^ 2 r, , , „ n 4 /■°° dP /i? r\ 

-P^-('-) = — -7^ \{r + a) pi^^ + 2 pisol + 5- / — i^ML — ,- di?, 

J r 

where the dimensionless kernel Kml is given in equation l l53t . 



(79) 



(80) 



2.4.8 Case of Diemand-Moore-Stadel anisotropy 

Finally for the anisotropy profile (eq. 1581 ) proposed bv lPiemand et alj j2004h . equation l l66b leads to (with Da from TableO 



GM , , 



2/3 



aVii r2 v;-Piso ^^^2)i/3 ^1/3 _ ^1/3 

for r < a, and to the radial solution (eq. 1741 ) for r > a. 



Piso + 



I/tt 



(ai/3 



(81) 



2. 4. 9 General form of the mass profile for specific anisotropy profiles 

Inserting equations ( 112b and J62b into equation J66b . one can obtain a general form for the mass profiles for the constant anisotropy, Osipkov- 
Merritt, Mamon-Lokas, and Diemand-Moore-Stadel anisotropy profiles: 

GM{r) _ 1 [°° \ Dp{r) 



-[3{r)]p(r)- 



VR^ 



dP 



R 



d^P 



dR - r2 



dR. 



(82) 



where, for the Diemand-Moore-Stadel anisotropy profile, the anisotropy radius a should be used for the upper integration limits. Equation l l82b 
allows to express the mass profile as a unique single integral of the observations, where Cp{r) and D/3{r) are given in Table |2l while the 
kernel Kjs is given by equations ( 137b . ( 148b . ( 153b . and ( 161b for the constant anisotropy, Osipkov-Merritt, Mamon-Lokas, and Diemand-Moore- 
Stadel anisotropy models, respectively. Equivalently, equation ( 182b can be used to formulate the circular velocity profile 
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vl(r) 



^ [1 - P{r)\ P{r) 



R 



Vi?2 - r2 dii2 



Cp{r)Kp{R,r) 



^AR2 



dP 
dP 



AR. 



(83) 



For isotropic models, equation J83t with = = if^ = recovers the second equality in equation J14b . 

In practice, writing the tracer density as p[r) = p(a) p{r/rs), where is the characteristic scale of the tracer, the projected pressure 
P{R/rs), equation l l83t yields 



as P{R) = P{r, 

-I 2 



crios(rs) 



S(r,)/[^r.p(r,)] 
[l-/3(r)] p(r/r,) 



X d^P Dp{r,x) 



VX2 



^dX2 



Ci3{rsx)Ki3 



dP 
dX 



dX . 



[nrsp{rs)] [1 - /3{r)] p{r/rs) 

fcosh"l(X„,ax/a:) 



O i(r x) ~ 
P"(x coshu) cosh II — - — \Cp(rsx) Kp(x coshu, x) sinhu — 1] P'(a:: cosh it) c??i (84) 

X 

where x — r/vs, X = R/vs, and where the second equality of equation ( 1841 ) is useful to avoid the singularity at X = a-, integrating out to 
the equivalent of, say, 10 r^, i.e. X^ax = 10 r^j/vs- All quantities on the right-hand side of the two equalities in equation iS4\ are known 
or assumed (the anisotropy profile). In particular, the numerator of the factor in front of the integral of the first equality of equation iS4\ is a 
function of the shape of the tracer density profile, found by Abel inversion (eq. (8)) of the surface density profile. 



3 TESTS 
3.1 Accuracy 

We test our mass inversion equations, on four anisotropy models: isotropic, constant, Osipkov-Merritt and Mamon-Lokas. For each o f 
these anisotropy models, we compute the projected pressure using equation ([6}, with the kernels given bv lMamon&LokasI ( I2005hll2006h . 
evaluated on a logarithmic grid from r — 0.01 to 100 in steps of 0.2 dex. The projected pressures P{R) were differentiated after 
cubic spline interpolation and the integral of equation J83l > was performed in steps of cosh^^(P/r) out to 100 Vs - We choose our mass an d 
anisotropy models by placing ourselves in the context of clusters of galaxies. We assume a one-component NFW model l lNavarro et aljl99^ . 
for which the dimensionless density and mass profiles can be expressed as 



pjxrs) _ (ln2- 1/2)' 
M{rs)/ {47vrl) ^ x{x + 1)^ 



—/ \ pv-^'s; 1^111 i- — J./ 



where is the scale radius, where the slope of the density profile is —2. We make no use of our assumption that the total and tracer density 
profiles are proportional. 

The anisotropy profile for dark matter particles in ACDM halos of the masses of clusters is close to the Mam on-Lokas model 
jMamon & L okas 2005b; Wojtak et al. 2008; Mamon et al. 200^ with anisotropy radius a ~ 0.18 r2oo dMamon & LokasI) or 0.275 r2oo = 
1.1 r a JlVlamon et alj) . We adopt a scaling of a = Ts for the Mamon-Lokas model and for the Osipkov-Merritt model as well, and we adopt 
a constant anisotropy model that is fairly radial but consistent with the anisotropy seen in ACDM halos: {3 = 0.4. 

Figure [T] shows the comparison of the circular velocity profiles obtained from the mass inversion equation ( I83t with the true circular 
velocity profiles. Despite the double differentiation of the projected pressure, the circular velocity (hence mass) profiles are recovered to a 
few percent relative accuracy or betterlfl except at the innermost point where the mass is overestimated by 4 to 10% in the four anisotropy 
models, because of the inaccurate cubic-spline interpolation of P(P) near the edges. The accuracy of the mass inversion is even better if we 
use a finer grid to measure the projected pressure before the cubic spline interpolation of P(P) and subsequent mass inversion: for example 
with the OM anisotropy, the maximum relative error in the recovered mass decreases with grid size from 9% (0.2 dex steps) to 0.6% (0.02 
dex steps). 



3.2 Robustness to small data samples 

We next test the accuracy of the recovered mass profiles when the data is sparse and noisy. We consider the case of velocity measurements 
in a cluster of galaxies. We assume that the cluster has 500 measured velocities within 5 (which is roughly the cluster virial radius), and 
assume for simplicity that we have line-of-sight velocity dispersions measured in 10 equal size radial bins centered from 0.25 to 4.75 r^. 
With A'^ = 50 velocities per bin, the velocity dispersions are known to a relative accuracy of ^l/2/(A-l) = 10.1% (e.g. lLupton|[l993h . 
and we fold this nois^ into the predicted line-of-sight velocity dispersion profile, using the same seed for the random number generator 
for all four anisotropy profiles. We extrapolate the projected pressures to larger radii by fitting a power-law to P(P) using the last 5 data 



° Both mass inversion and deprojection appear to be unstable at radii r < Vs for the Mamon-Lokas anisotropy model when a is exactly set to rs, when using 
our Mathematica routines (but this odd behavior is not present when tested with other software). The figure shows the case a = 1.001 r^. 
^ We neglect the noise on the uncertain surface density profile, which contributes negligibly to the noise in the projected pressure in comparison to the noise 
in the velocity dispersion. 
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Figure 1. Left: Adopted anisotropy profiles: isotropic (solid black line), Osipkov-Memtt (eq. 1451 . witli a = rs, dotted red line), Mamon-Lokas (eq. 1491 . 
with a = Ts, short dashed green line), and /3 = est = 0.4 {long dashed blue line). Right: Accuracy of the mass inversion (in the absence of noise): ratio of 
inferred (eq. 1831 . using eq. (6) to first evaluate P on logarithmic grid of 0.2 dex steps, and using the dimensionless functions of Table[2] and the dimensionless 
kernels of equations 1481 . 1531 . and 1371 . for the latter three anisotropy models) over true NFW mass profiles for the four anisotropy models shown in the left 
panel. 

points, at outer linearly spaced outer radii, with the same spacing as the data, and then fit a 4th order polynomial to the set of observed and 
mock-extrapolated data. We repeated these tests 5 times with different seeds for the random generator. 

The left-hand panel of Figure |2] shows the accuracy of the mass inversion is much worse than in the academic case with no noise. In 
particular, the extrapolation errors at radii lower than the lowest data point make the inner mass profile inaccurate to factors often greater 
than 2. With the isotropic, (3 = 0.4, and Mamon-Lokas anisotropy models, the mass profile is nevertheless recovered to typically better than 
20% accuracy for r > 0.8 Ts, out to twice the radius of the last data point. However, the large error bars show that there is a large scatter 
in the accuracy of the recovered mass profile for different randomly generated projected pressure profiles. The right-hand panel of Figure|2] 
gives the second highest error among the five tests performed, for each given radius and anisotropy model. Typical such 80-percentile errors 
are of the order of 70% for r > 0.8 r^. Surprisingly, this typical error decreases to only 20% at high radii (r > 8 r^), despite the fact that the 
projected pressure is extrapolated beyond r — 4.5 r^. 

3.3 Robustness to the wrong anisotropy model 

The essential ingredient to the mass inversion is the knowledge of the velocity anisotropy profile. How wrong can the mass inversion be if 
the incorrect a nisotropy profile is u sed? We adopt the Mamon-Lokas anisotropy model with a = Vs similar to what is found for cluster-mass 
ACDM halos jMamon et al]|2009l) and compute the projected pressure for an NFW model with this anisotropy model. We then perform the 
mass inversion assuming other anisotropy profiles to see how off we are. In this exercise, we assume perfect data, i.e. no noise. 

The left panel of Figure[3]shows that the mass profile is recovered to better than 33% accuracy for all anisotropy models at r > 4 rs, i.e. 
beyond the virial radius. Within the virial radius, the Osipkov-Merritt underestimates the mass by as much as a factor 3 around 2rs, but is 
much more precise at very low radii. The P = 0.4 model is accurate for r > r^, as expected as it resembles there the Mamon-Lokas model, 
but underestimates the mass by increasingly large factors at radii r < Vs, and the recovered mass actually goes negative at r < 0.17 r^. 
Finally, the isotropic model finds the correct mass to within 30% at all radii, usually overestimating the true mass. Interestingly, at r ~ 7rs, 
all four anisotropy models lead to the correct mass to within 5%. 

The right panel of Figure [3] indicates that the recovered mass is not very sensitive to the assumed anisotropy radius, as the mass is 
recovered to 20% accuracy, unless the anisotropy radius is assumed to be 10 times lower than it actually is. This graph also shows that 
at r ~ 3rs (i.e. roughly two-thirds the virial radius of clusters), the mass is correctly recovered to better than 5% for our five choices of 
anisotropy radius. 



4 DISCUSSION 

The mass inversion algorithm presented in this work generally (eq. 1691 ), requires two steps; 1) deprojection and 2) inserting the radial 
pressure in the Jeans equation to derive the mass. The deprojection (eq. 1321 ) requires a single integral involving of a quantity, piso (eq. |10| ) 
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Figure 2. Robustness of the mass inversion to small data samples. Left: Same as right panel of Fig.[T] but for the projected pressure profile measured on a 
linear grid of 10 radial bins from 0.25 to 4.75 rs, with 20% relative gaussian errors on the projected pressure (i.e. 10% errors on the line-of-sight velocity 
dispersion measurements based upon 50 velocities per radial bin). The error bars show the standard deviations on 5 tests with different seeds for the random 
number generator. The points and eiTor bars are slightly shifted along the x-axis for clarity. Right: 2nd highest error out of 5 tests on recovered mass profile. 
A value of unity indicates a perfect recovery of the mass. 



that is itself a single integral involving the derivative of the observed projected pressure. The second step (mass inversion) also requires a 
single integral involving the derivative of piso- Each differentiation of the data introduces errors, and the full mass inversion requires three 
single integrals. For the special cases of simple anisotropy models, we find it preferable to write the mass profile with a single integral 
involving the double derivative of the observed projected pressure. Indeed, this requires a single smoothing operation before differentiation, 
thus leading to more accurate results, even if the mathematical formulation of the deprojections and mass inversions for each of the simple 
anisotropy models has stro ngly increased the number of equations in this article. 



While this work (andlWolf et alj (|2009|) ) used the Abe l inversion for th e kinematic deprojection, one can alternatively apply Fourier 



methods (see also lKalal & Nugentll 19881 and Kalnais cited in ISahaetalJll996h . Indeed, structural and kinematic projection can be written as 
a convolution: 

/oc 
f{x)K(X~x)Ax, 
-oo 

where X = E?,x = r^, F(X) and f{x) co rrespond to either and p(r) (structural projection) or P{R) and (1 — /?) p + (]pds/s 
(anisotropic kinematic projection, following IWolf et al.L which simplifies to p for isotropic kinematic projection), and where 

\ y > . 

Hence, with the convolution theorem, deprojection is obtained by applying an inverse Fourier transform to 

fiLo) = ^ ^ = (1 + ^) sgn(^) 7f F(a;) ^ [1 + ^ sgn(..)] , 
V27r K{iu) V 27r y Itv 

where /(a;), K{uj) and F{uj) are the Fourier transforms of f{x), K(y) and F {\X\), respectively (note the absolute values in the last term). 
A comparison of the accuracy of the two deprojection techniques is beyond the scope of the present article. 

Our mass inversion algorithm should serve as a useful technique to get around the mass-anisotropy degeneracy in the case where the 
anisotropy profile is thought to be known. As mentioned in the end of Sect. I, there is a good convergence on the anisotropy profiles of ACDM 
halos as well as those of elliptical galaxies formed by binary merge rs of spiral galaxies. Mo reover, the anisotropy profile in many simulations 



appears linearly related to the slope of the (tracer) density pr ofile i Hansen & Moordl 20061) . and this can be used to lift the mass-anisotropy 
degeneracy. A first application of our algorithm was given bv lBiviano & Saluccil ^200^ for the analysis of stacked clusters of gala 



ilaxies. 
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Figure 3. Robustness of the mass inversion to tlie wrong clioice of anisotropy protile. Left: Same as riglit panel of Fig. [T] where the true anisotropy profile 
is now always an a = rs Mamon-Lokas model, but assuming that of the other three anisotropy models. Right: Same as left panel, but where the assumed 
anisotropy profile is always the a = rs Mamon-Lokas model, but with 5 different choices for a/r^. 

The ma.ss inversion technique ha.s the advantage of producing a non-parametric0mass profile, which can then be used to test the popular 
parametrizations of the mass profile (or alternatively of the density profile, the circular velocity profile or the density-slope profile). 

In Sect. [3] we show that, for a mock NFW galaxy cluster with mildly increasing radial velocity anisotropy as seen in ACDM halos 
and with typical line-of-sight velocity dispersion profiles, measured with 50 velocities per radial bin, the mass inversion should be accurate 
to typically better than 70% relative errors at most radii and better than 20% for anisotropy models other than the Osipkov-Merritt one at 
r > Sts. The relatively high errors are a consequence of the double derivative of the observed projected pressure, d?P/AB?, that enters 
the mass inversion equation ([68} or ([69}, through the term piso, or in equation ( 1821 ) or ([83} for the special cases of anisotropy models. The 
errors are high at radii smaller than the first radial bin of the observed line-of-sight velocity dispersion profile. This illustrates the concept 
that kinematical modelling can only recover the mass and anisotropy at radii corresponding to the projected radii of the data. Nevertheless, 
with power-law extrapolations of the data to outer radii, we show that the mass inversion can recover mass profiles with good accuracy far 
beyond the outermost data point. Note that the mass inversion involves integrals out to infinity (e.g. eq. 1831 ). so one expects that the method 
should be most accurate when the tracer density profile falls fast at large radii. Our use of the NFW model for the tracer, with its shallow 
outer slope of d In p/dln r = —3 is thus expected to provide poorer results for the mass inversion than for steeper tracer density profiles. 

We found that the recovered mass is correctly returned, independently of the shape of the anisotropy profile at r = 7rs, and indepen- 
dently of the anisotropy radius for our chosen a nisotropy model a t r — 3 rs. A similar independence of the recovered mass on the assumed 
anisotropy profile has been recently noticed by Wolf et al, (2009h in the context of dwarf spheroida l and ellipt ical galaxies (for which the 
dark matter may not follow the stars, which themselves do not follow the NFW model). However, IWolf et alj prove analytically that this 
robustness to the ani sotropy mod el occurs near the the radius of slope —3. Now, the NFW model has shallower slopes everywhere, reaching 
—3 at infinite radius. IWolf et al J notice that, for density profiles similar to those of ellipticals and dwarf spheroidals, the radius of slope —3 is 
close to the half-mass radius. In contrast, in the current context of clusters, the NFW model is divergent in mass (eq. 1861 ), and the concept 
of half-mass radius is ill-defined. Moreover, the radius where the mass is recovered for all anisotropy models tested is at 7 scale radii, which 
is outside the virial radius, hence not comparable to the h alf-light radius of elliptical and dwarf spheroidal galaxies. Fixing the anisotropy to 
the Mamon-Lokas model (which lMamon & Lokasl2005bl found to be a good fit to the anisotropy profile of the halos in ACDM cosmological 
simulations), the recovered mass is most robust to the anisotropy radius at 3 Vs, which is roughly two-thirds of the cluster virial radius, again 
not directly comparable to the half-light radius of dwarf spheroidals and ellipticals. 

The mass inversion technique is thus a useful complement to the set of tools one has to lift the mass-anisotropy degeneracy in spherical 
systems. Mass inversion is certainly not the privileged tool, but should be considered as one of many tools for the exploratory data analysis of 
spherical systems viewed in projection, in addition to anisotropy inversion, fitting models to the line-of-sight velocity dispersion and possibly 
kurtosis profiles, and fitting models, distribution functions, orbits and A'^-body systems to the distribution of particles in projected phase 
space. Ideally, one would analyze the kinematics of spherical systems using a variety of these tools. We are preparing such global analyses 
on dwarf spheroidal and elliptical galaxies, as well as on groups and clusters of galaxies. 

^'^ Although the algorithm uses a parametrized anisotropy profile and a smooth fit through the projected pressure profile, the mass profile that comes out is 
non-parametric. 
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APPENDIX A: ABEL DEPROJECTION 

In this appendix, we remind the reader of the derivation of the deprojection of equation Q with the Abel inversion. Consider 
Replacing E(i?) in equation JAlt by its definition in equation jT), one finds, after inverting the order of integration: 

J(r) = 2 / p{s)sAs [ T7?- (A2) 

The internal integral in equation l lA2t is equal to tv/2, as inferred from the substitution sin'^ 6 — {R^ — r^)/ (s^ — r^). Hence, 

/> oo 

J(r) = TV p(s) sds , 

J r 

and therefore 

Pir) = ~-^. (A3) 
nr ar 

Now integrating equation l lAlb by parts, one gets 

dE / 2 2\l/2 



J(r)= lim y/R2 - r2 E(i?) - / ^ (r' - rY dR , (A4) 

J r 

For all realistic density profiles, falls faster than R~^, as is the case for the surface density profiles of globular clusters, elliptical 

galaxies and clusters of galaxies. Hence, the surface term in equation ( IA4b is zero and one can write 

dj^ r dE dR 

dr X dR (i?2_ ^2)1/2 • ^ ^> 

Inserting the derivative of J of equation l lA5l l into equation JA3b leads to equation {Sjl. The surface term that survived when E oc 1/R 
disappears in the derivative. 



APPENDIX B: KINEMATIC DEPROJECTION FOR THE TANGENTL4L DYNAMICAL PRESSURE 

In this appendix, we derive equations l l31t and ( 132b for the tangential dynamical pressure. 
Differentiating equation ( 123b . one finds to get the differential equation 

I P pe I 

Pe--, s — =Piso. (Bl) 

I — p r 

Now, if we write 

, /3 pg _ 1 djfpe) 

Pe — ^ — — T — A I 

1 — p r / dr 

then equations dBlb and ( IB2b lead to 

Pe{r) = -jr-TT / fpLods, (B3) 

where the upper limit at infinity ensures that pg = {1 — (3) pa^ does not reach negative values at a finite radial distance. But equation ( IB 2b 
directly gives 
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din/ _ I3{r) 
dlnr ^ 1 - /3(r) ' 

hence 

g{r) = 5(ri) exp [ -] (B4) 



1-/3 s 

for any arbitrary ri. With equation JB4b , equation JB3I > allows one to recover equation Olt : 
pe{r) = -exp / ^— / exp - / ^ — Piso ds , 



X , 1 - /3 s / / V / 1 - /3 i , 

cxp(^-£^^^ pLod., (B5) 

where the second equality is obtained adopting r\ — r. 

Integrating by parts the integral in equation JB5b . we finally recover equation l l32b 

, , r p ( r p At\d8 
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